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Abstract 

A paradigm for isothermal, mechanical rectification of stochastic fluctuations is in- 
troduced in this paper. The central idea is to transform energy injected by random 
perturbations into rigid-body rotational kinetic energy. The prototype considered in 
this paper is a mechanical system consisting of a set of rigid bodies in interaction 
through magnetic fields. The system is stochastically forced by white noise and dis- 
sipative through mechanical friction. The Gibbs-Boltzmann distribution at a specific 
temperature defines the unique invariant measure under the flow of this stochastic 
process and allows us to define "the temperature" of the system. This measure is also 
ergodic and strongly mixing. Although the system does not exhibit global directed mo- 
tion, it is shown that global ballistic motion is possible (the mean-squared displacement 
grows like t^). More precisely, although work cannot be extracted from thermal energy 
by the second law of thermodynamics, it is shown that ballistic transport from thermal 
energy is possible. In particular, the dynamics is characterized by a meta-stable state in 
which the system exhibits directed motion over random time scales. This phenomenon 
is caused by interaction of three attributes of the system: a non flat (yet bounded) 
potential energy landscape, a rigid body effect (coupling translational momentum and 
angular momentum through friction) and the degeneracy of the noise/friction tensor 
on the momentums (the fact that noise is not applied to all degrees of freedom). 

1 Introduction 

Mechanical rectification was introduced to describe devices that convert small-amphtude 
mechanical vibrations into directed rotary or rectilinear mechanical motion [6]. The pur- 
pose of this paper is to extend this concept to mechanical systems in isothermal environments 
subjected to stochastic fluctuations. Although the second law of thermodynamics prevents 
directed mechanical motion from thermal fluctuations in an isothermal environment, vio- 
lations of the second law of thermodynamics on certain time-scales in microscopic systems 
have been theoretically predicted [10; 12] and experimentally observed [7; 20]. 

In this paper we exhibit a system characterized by ballistic non-directed motion at 
uniform temperature. Set /i to be the unique invariant measure of the system and x(t) the 
displacement of the system at time t then ^ a.s. limt^oo{x{t) — x{0))/t but E^[{a;(i) — 
E(x(i))}^] ^ P. Moreover the system is characterized by "dynamic" meta-stable states in 
which it experiences ballistic directed motion over random time scales (flights) in addition 
to "static" meta-stable states where the system is "captured" in potential wells. 

The mechanism introduced in this paper is different from the one associated to the 
Gallavotti-Cohen fluctuation theorem [10]. It is based on the fact that one can obtain 
anomalous diffusion by introducing degenerate noise and friction in the momentums. This 
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anomalous diffusion manifests itself in the just mentioned flights over time-scales which are 
random and the probability distributions of the flight-durations have a heavy tail. Moreover 
we expect that it will be of relevance to Brownian motors [2; f 9] since it shows how thermal 
noise can be used to get ballistic transport, in other words how to obtain fast and efficient 
intracellular transport without using chemical energy. Although global ballistic transport 
is achieved it is not directed and hence work is not produced from thermal fluctuations 
(this mechanism is not in violation of the second law of thermodynamics). However this 
mechanism would be sufficient for intracellular transport of a large fraction of the material 
being carried to target areas with ballistic speed without draining cellular energy reserves. 




Figure f .f : Picture of mechanical system. The physical system consists of a magnetized axisymmetric 
top on a flat surface and a magnetized ring as shown above. The dipole moments in the ring are oriented 
radially. 

The prototype analyzed in this paper consists of a magnetized top (ball) on a surface 
interacting with a suspended magnetized ring as shown in Fig. f .f . The following dynamics 
is observed: when one lowers the magnetic ring to within a certain range of heights and 
then tilts the ring, one observes the top transition from a state of no spin about its axis of 
symmetry to a state of nonzero spin. The reader is referred to the following urls for movies 
and simulations of the phenomenon: 

http : //www. acm. caltech. edu/~nawaf /BallisticTransport/ 
http : //www. acm. caltech. edu/~owhadi/BallisticTransport/ 

This phenomenon seems counterintuitive and non-Hamiltonian, since one would expect the 
angular momentum of the top to be conserved (since without friction the Hamiltonian of 
the system is invariant under an rotation of the ball). The origin of this mechanical 
device can be traced back to the work of David Hamel on magnetic motors in the unofficial 
sub-scene of physics [f5]. The mechanism presented in this paper when the device is not 
at uniform temperature could in principle be used as a method of extracting energy from 
macroscopic fluctuations. 

2 Preview of the Paper 

In section 4 Hamel's device is analyzed through an idealized model based on magnetostatics 
and the spinning of the top is caused by the introduction of surface frictional forces. Sim- 
ulations (figure 1.2) are done using variational integrators and concur with experimental 
observations. 

In section 5 the magnetic ring is allowed to be dynamic, and a fixed outer ring of a finite 
number of magnetic dipoles is introduced to stabilize it (see figure 5.1). The inner magnetic 
ring is excited through white noise applied as a torque. The steel ball and inner ring are 
coupled through a magnetostatic potential and the motion of the steel ball is dissipative 
through slip friction. 
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Figure 1.2: Snapshots of Nonconservative Tilted Simulation. The above are snapshots of the 
simulation described in Fig. 4.9. The top is initially set with its axis of symmetry pointing nearly vertical. 
The axis of symmetry then aligns with the ambient magnetic field. However, the state is not a minimum of 
Ve, and hence, the state is unstable. The top moves towards a state that minimizes its magnetic potential 
energy, and acquires spin in this process. 



The mathematical description of the fluctuation driven motor is obtained by generahzing 
Langevin processes from a system of particles on a linear configuration space to rigidified 
particles whose configuration space is the Lie group SE(3). The question of random pertur- 
bations of a rigid body was treated by previous investigators who added perturbations to 
the Lie-Poisson equations without potential or dissipative torques [13; 14]. Hence, they do 
not consider generalizing Langevin processes to the Lie-Poisson setting. 

We refer to figure 2.1 for a plot of the angular position of the ball versus time for different 
values of noise amplitude a. The simulation of the system exhibits two distinct kinds of 
metastable states. In the first kind the ball is "stuck" in a magnetic potential well whose 
depth depends on the position of the inner ring, and moves into another potential well when 
the energy barrier between the two wells is close to minimum (a phenomenon known as 
stochastic resonance) or transitions to a metastable state of the second kind. In the second 
kind the ball spins in circles clockwise or counter-clockwise in a directed way for a random 
amount of time (that numerically and heuristically observe to be exponential in law) until 
it gets stuck in a potential well. 




(a) a = 0.0002 (b) a = 0.0002 (c) a = 0.00025 



Figure 2.1: Angular-position of magnetic ball (non-uniform temperature). The angular com- 
ponent of the center of mass is plotted for three different realizations. The plots show the magnetic ball 
transition between two meta-stable states: noise-driven and inertia-driven motion. 

Since the system associated to figure 2.1 is described by a generalized Langevin process 
it is possible to introduce a notion of temperature defined as ("°'^° amplitude) ^ However, 

^ ^ rriction constant ' 

since noise but no friction is applied at the level of the inner ring and friction but no noise 
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is applied at the level of the ball, the system is not at uniform temperature, i.e., the inner 
ring is at infinite temperature and the ball is at zero temperature. Nevertheless, one could 
in principle use such a mechanism as a way to extract energy from macroscopic fluctuations. 

In a second step frictional torque is introduced to the inner ring and thermal torque 
(white noise) to the ball, so that the generator of the process is characterized by a unique 
Gibbs-Boltzmann invariant distribution. Throughout this paper we interpret this property 
as placing the mechanical system at uniform temperature. We refer to figure 2.2 for a plot 
of the angular position of the ball versus time for different values of noise amplitude a (and 
temperature). The system is still characterized by ballistic directed motion meta-stable 
states. 




(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 



Figure 2.2: Angular-position of magnetic ball (uniform temperature). Four different realizations 
of the angular component of the center of mass are plotted. The plots for a = 0.0002 show the magnetic 
ball transition between two metastable states: noise-driven and inertia-driven motion. If a is increased to 
a = 0.00025, realization (c) shows a similar transition. Realization (d), however, does not show such a 
transition, i.e., the ball never transitions to a metastable state of the second kind. 




Figure 2.3: Sliding Disk. Consider a sliding disk of radius r that is free to translate and rotate on a 
surface. We assume the disk is in sliding frictional contact with the surface. The configuration space of the 
system is SE(2), but with the surface constraint the configuration space is just M X SO(2). 

To understand the behavior of the fluctuation driven motor prototype, the paper con- 
siders in section 3 a sliding disk (figure 2.3) that has the same essential behavior, but whose 
configuration space is SE(2). The solution of this simplified system is a Langevin process 
with degenerate noise and friction matrices in the momentums. The disk is free to slide and 
rotate. Assume that one rescales position by r and time by some characteristic frequency 
of rotation or other time-scale. The dimensionless Lagrangian is given by the difference of 
kinetic and potential energy 



L{x,v,9,iu)^-v^ + -u;^-U{x) 
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U non-constant => Ballistic Diffusion 




U constant Nornnal Diffusion 



Figure 2.4: Ballistic vs. Normal Diffusion. If U is non-constant theorem 3.1 implies tliat tlie mean 
squared displacement with respect to the invariant law is ballistic. Numerically, flights are observed in 
the x-displacement as shown in the diagram. However, if U is constant this diffusion is normal, i.e., the 
i-displacement behaves like Brownian motion. 



where [/ : M ^ M is assumed to be smooth and periodic. \i U = cos(a;), figure 2.5 shows 
that the shding disk is a modified one-dimensional pendulum. The contact with the surface 
is modelled using a sliding friction law. For this purpose we introduce a symmetric matrix 
C defined as, 

^ " [l/a 1/^2 ■ 

C is degenerate since the frictional force is actually applied to only a single degree of freedom, 
and hence, one of its eigenvalues is zero. In addition to friction the system is excited by 
white noise so that the governing equations become 



dx 
dO 



vdt 
Lodt 



dv 








V 




dBy 




dt-cC 






doj 













dB^_ 



(2.1) 



where C^/^ is the matrix square root of C. Let E denote the energy of the mechanical 
system given by, 

U{x), (2.2) 



1 9 1 9 



and let f3 = 2c/q;^ 

In [4], we show that the Gibbs-Boltzmann distribution 

H = exp {-PE) , 

defines the unique invariant measure under the fiow of the sliding disk. This measure is also 
ergodic and strongly mixing. Using this result we prove in § 3, that if U is non-constant then 
the x-displacement of the sliding disk is fi a.s. not ballistic (cf. proposition 3.2). However, the 
mean-squared displacement with respect to the invariant law is ballistic (cf. theorem 3.1). 
More precisely, we show that the squared standard deviation of the x-displacement with 
respect to its noise-average grows like . This implies that the process exhibits not only 
ballistic transport but also ballistic diffusion. If U is constant then the squared standard 
deviation of the cc-displacement is diffusive (grows like t) . 

In the numerics the sliding disk is initially at rest and averages are computed with 
respect to realizations. Numerically one observes the following consequences of this ballistic 
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behavior. Figure 2.6 shows that when U is non-constant then the motion is characterized 
by meta-stable directed motion states. Figure 2.7 shows that the mean square displacement 
E[a::(t)^] grows hke t (with respect to time) when U is constant or in the case of a one 
dimensional standard Langevin process (control) whereas it grows like as soon as U is 
non-constant (the motion becomes ballistic). A diagram comparing the solution behavior 
in the U constant and non-constant cases is provided in Figure 2.4. Figure 2.8 corresponds 
to the plot of 

Cov(x(t + 2s) - x{t + s),x{t + s) - x{t)) 
{Va.T{x{t + 2s) - x{t + s)))HVar{x{t + s) - x{t)))h 

for t large as a function of s. It clearly shows that the system is characterized by long time 
memory/correlation when U is non-constant whereas it has almost no memory when U is 
constant or in in the case of a one dimensional standard Langevin process. 




Figure 2.5: Ballistic Pendulum. If the dimensionless potential is U = cos(a;), then the sUding disk is 
simply a pendulum in which the bob in the pendulum is replaced by a disk and the pendulum is placed 
within a cylinder as shown. 
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(a) U-symmetric (b) U-asymmetric (c) U-flat (d) Control 

Figure 2.6: Sliding Disk at Uniform Temperature, h = 0.01, a = 5.0, c = 0.1. The mean of the 

x-displaccmcnt of the disk for U symmetric, asymmetric, flat and control case. Figures (a) and (b) show that 
when U is non-constant, directed motion as a incta-stable state is jjossible. On the other hand, (c) and (d) 
do not show such behavior. The figure on the top superposes these graphs in a single plot for comparison. 




10' 10' 10* 10' 10' 



Figure 2.7: Sliding Disk at Uniform temperature, h = 0.01, a = 5.0, c = 0.1. A log-log plot of the 
mean squared displacement of the ball. It clearly shows that the x-position exhibits anomalous diffusion 
when U is symmetric or asymmetric. In the control and flat U cases the diffusion is normal. 
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(a) U-symmetric (b) U-asymmetric (c) U-flat (d) Control 

Figure 2.8: Sliding Disk at Uniform Temperature, h = 0.01, a = 5.0, c = 0.1. The correlation of the 
X displacement of the disk for U symmetric, asymmetric, flat, and the control case. When U is symmetric 
or asymmetric the correlation in the x displacement is nonzero for a certain time-scale which is larger than 
the characteristic time-scales associated with U and with the friction factor. However, in the other cases 
the correlation is negligible. The figure on the top superposes these graphs in a single plot for comparison. 



3 Sliding Disk: Simplified Model 

To understand the behavior of the prototype stochastic mechanical rectifier which will be 
treated in detail in subsequent sections, we designed a simplified model whose configuration 
space is SE(2). The system consists of a disk sliding on a surface as shown in Figure 2.3. 
The effect of the outer ring is modelled as a periodic, one-dimensional potential; while the 
effect of the inner ring is incorporated into white noise. The SDK for the isothermal, sliding 
disk is a Langevin process, but the noise and friction matrices are degenerate. A statistical 
numerical analysis discussed below shows that this process has two interesting statistical 
properties which are both linked to rigid-body inertial effects: 1) a non-trivial correlation 
on certain time-scales and 2) ballistic diffusion. 

Ballistic Transport at Uniform Temperature The proof that the sliding disk is at 
uniform temperature is based on finding the generator for (2.1) and showing that the Gibbs 
measure is the unique, invariant measure under the flow of this generator. 

Theorem 3.1. Let x:=TrxMxTxM denote the phase space of the sliding disk (T standing 
for the torus of dimension one). Set :— {x,v,9,uj) to he the solution of (2.1). Let E denote 
the energy of the sliding disk given by, E = \v'^ + \<yLu'^ + U {x) . Set (3 = 2c/a^ and let ji he 
the Gibbs probability measure defined by 

M(dO := -^dC (3.1) 

where Z -.^ J e'^^^d^. Lf U is non-constant, then the Gibbs measure fj, is ergodic and 
strongly mixing with respect to the stochastic process ^. Furthermore, it is the unique, 
invariant probability measure for the stochastic process ^. 

This result is a special case of the proof provided in [4]. To determine the measure is 
invariant, one computes the infinitesimal generator C associated to the stochastic process ^ 
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and proves that if f is fi - measureable then, 

/ Cf^i{d0 = 



Using the fact that the measure is ergodic, one can readily prove that the x-displacement 
itself is not ballistic. 

Proposition 3.2. Provided that U is non-constant, then fi a.s. 

lin. '-i^hj^ ^ 0. 

t — >oo t 

Proof. Since is ergodic with respect to ^, 



t^oc t 



lini - / vdt^ I v^i{d^) = 0. 
t-.oo t Jq 



Moreover, as stated in the following proposition, the squared standard deviation of the 
X + ^-degree of freedom of the sliding disk grows like t. Let E denote the expectation with 
respect to the Brownian noise and the expectation with respect to the Brownian noise 
and the initial configuration (sampled from the invariant measure /i). Set Xt — x{t) and 
0t = 0{t). 

Proposition 3.3. The squared standard deviation of the Xt + 9t-degree of freedom is diffu- 
sive, i.e., 

E^[ixt + et~E[xt + et])^] 2«V2 



lim ^ 

t — >OQ 



t 



C2(cr2 + 1) 



(3.2) 



Proof. First, diagonalize the diffusion and friction matrices in (2.1) using the following 
invertible matrix: 



1 



l + a 



-1 a 
1 1 



as follows 
V 



dv 
duj 



=V 







dt — c 





^ 



V 



dt 



a{a + l) 



Simplifying this expression yields the following pair of equations: 

d {—V + auj) = dxUdt, 

d{v -\- uj) = —dxUdt — C7 (w + w) + a{dBy + dB^). 



"0 


0' 


V 


'dB,' 





1 


dB^_ 



(3.3) 
(3.4) 



where 



(cr + 1) _ a{a+l) 

7 = > Oi ■ 



Set Bs := {By + B^)/V2. Integrating (3.3) and (3.4) gives 

—v-\-auj — —vo + o'uJo+ / dxUds, 

Jo 

v + iJ = {va+ uJo)e~''^' - [ e-"^^'-"^ d^Uds + V2a [ e-^^^'-'^dBs 

Jo Jo 



(3.5) 
(3.6) 
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Integrating (3.6) gives 

{xt + Ot) - (xo + 6*0) 



C7 Jo C7 

f* 1 _ g-c7(t-s) 



dxUds 



C7 



The result follows by observing that 



dxUds = {-vt + CTWt) - (-Wo + cwo) 



(3.7) 



However, one can prove that the squared standard deviation of the —Xt + aOt-degiee of 
freedom grows like . This implies that the process exhibits not only ballistic transport but 
also ballistic diffusion along this degree of freedom. 



Proposition 3.4. Assume that U is non constant, then 



E,, 



lim sup 

t — >oc 



Xt + a9t - V\-xt + aOt])' 



t2 



< 4 



1 + g 



and 



lim inf ■ 

t — >oo 



xt + cjet-¥.[-xt + (Tet]y 



> 



1 l + cr 



t^ - A 13 

Proof. From Cauchy-Schwartz inequality one obtains that 

— + (Tcij,) ds\ < — 

Hence 

E„ \( -xt + (jet - (-.Tn + (76'n))^l 

< A* [(-Wo + CTCJo)^ 



^ j {-Vs + cFLOs)ds^ <-j: [ {~Vs + cruis)'^ ds. 



lim sup 

t — >oc 



t^ 



We obtain the first inequality of the proposition by observing that 

l + cr 



m[(-wo + o-t^o)^] = 



and 



E„ 



<E„ 



{E[-xt + aOt] - i-xo + a9o)y 
Let us now prove the lower bound. Integrating equation (3.5) gives 

/■* s 

-xt + aOt - E[-a;t + aOt] ^t (1 - -)(a,C/ - E[d,U])ds 

Jo t 



Write 



At= f {l-^){dxU-E[dxU])ds 
Jo ^ 



(3.8) 



(3.9) 



(3.10) 



{-xt + Get-{-xo + aea)f (3.11) 



(3.12) 



(3.13) 
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and ^ 

Bt= f ^{d^U -nd^U])ds (3.14) 
Jo 

Observe that 

At + Bt^ i-vt + aujt) ~ E[-vt + aut]. (3.15) 
Since /i is strongly mixing when U is non constant, it follows that [18] 

lim inf E^iiAt + B^f] = fi[i-vt + Guj,f\ (3.16) 

Furthermore 

E^[(A + B,f\ < 2(E^[A?] + E^[Bf]) (3.17) 

and since the law of the process {xt, 0t,Vt,i-Ut) remains invariant under by reversing time 
and flipping the velocities Vt,LUt we deduce that Ep[yl(] = Ep[i3('^] and 

lim inf E^[A2] = l^i[{-vt + acJtf] (3.18) 

We conclude by the taking the expectation of square of (3.12) with respect to E^. ■ 

The following theorem is a straightforward consequence of the previous propositions. 
Theorem 3.1. We have 
• If U is constant then 

E4(x,-EH)^] 2a^a^ 

hm = „, „ — --r (3.19) 



• If U is non constant then 

r E4(x,-EM)^] 4 

hm sup < 3.20 

t^oo t-^ ^(1 + ct) 

and 

lim inf ^^[(-^-n^^^n > / ^ (3.21) 
Remark 3.1. Using Ito's formula and (2.2) one obtains that, 

dE = —c{v + uj^dt + — ( 1 H — ) dt + martingales. 



Integrating this expression gives 

E{t) - E{Q) = -c J iv + (^Yds + Y + ~^ ^ + martingales. 

The first and second terms represent the energy loss due to friction and the energy injected 
due to the noise respectively. 

Remark 3.2. Setting vt — v{t), it follows from (3.5) and (3.6) that 

vt= — / — d^Uds 

+ 1 Jf) cr + 1 



V2a 



'^-f(t-')dB,. (3.22) 



The long term memory effect exhibited in Figure 2.8 has its origin in the term dxU{x{s))ds 
in (3.22). That term is equal to zero when U is constant, and itself has its origin in the 
rigid body interaction between rotation and translation and the fact the friction matrix is 
singular. 
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Stochastic Variational Integrator To simulate the dynamics of the sliding disk at 
constant temperature, a stochastic variational Euler method is applied [5]. The discrete 
scheme for the isothermal case is given explicitly by: 



Xn+l = Xn + hVn+1, 
(^n+l — Sn + huJn+1, 



Vn+l 




Vn 


+ h 


-dxU{Xn) 


-hcC 


Vn 










^n 









cruon 







(3.23) 



It is an explicit, first-order strongly convergent method. 

Simulation We will consider four different systems to simulate. The first three are sliding 
disks with a symmetric, asymmetric, and flat potentials: 



U{x) 




sm{x) 


symmetric 


sin(x) + 0.4sin(2a;) 


asymmetric 





flat 



The fourth case is a control consisting of a simulation of a 1-D Langevin process at the same 
temperature and with a symmetric potential: 

dX = Vdt, 

dV = -cVdt - cos{X)dt + adBv- 

For all of the simulations the sliding disk is initially at rest, r = 0.25, m = 1.0, J — mr^/2 
and the friction and noise factors are as indicated in the figures. The means of the x 
displacement of the disk as shown in Fig. 2.6 are very small compared with the spread as, 
e.g., shown in the histogram of the final position of the ball as shown in Fig. 3.3. However, 
the mean squared displacement shows ballistic diffusion in the cases when U is symmetric or 
asymmetric and normal diffusion otherwise (see Figures 3.1-2.7). The time-scale associated 
with this ballistic diffusion is plotted in Fig. 2.8 which shows the correlation in the x- 
displacement when U is symmetric or asymmetric. This time-scale is much greater than 
the characteristic time-scale associated with the friction or the potential. Recall from the 
integral expression of the velocity, that when U is zero a rigid-body term is neglected. This 
demonstrates the important role of the rigid-body effect in the ballistic diffusion of the 
^-displacement when U is symmetric and asymmetric. 

Finally we also consider adding a non-degenerate, but anisotropic dissipation matrix to 
the sliding disk. It is numerically observed that if the anisotropy is large enough, E(a;i) is 
ballistic for a long period of time. 
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(a) U-symmetric (b) U-asymmetric (c) U-flat (d) Control 

Figure 3.1: Sliding Disk at Uniform Temperature, h = 0.01, a = 5.0, c = 0.1. Prom left: the 
mean squared position of the disk for U symmetric, asymmetric, flat, and the control. The figure on the top 
superposes these graphs in a single plot for comparison. 




(a) U-symmetric (b) U-asymmetric (c) U-flat (d) Control 



Figure 3.2: Sliding Disk at Uniform Temperature, h = 0.01, a = 5.0, c = 0.1. From left: "diffusion" 
of x-displaccmcnt of the disk for U symmetric, asymmetric, flat, and the control. For the cases when U is 
flat and the control, the diffusion is normal. Whereas in the other cases the diffusion is ballistic. The figure 
on the top superposes these graphs in a single plot for comparison. 
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(a) U-symmetric (b) U-asymmetric (c) U-flat (d) Control 

Figure 3.3: Sliding Dislc at Uniform Temperature, h = 0.01, a = 5.0, c = 0.1. From left: histogram 
of the x-displacement of the disk at T = 50000 for U symmetric, asymmetric, and flat. We observe a wider 
spread in the cases when U is symmetric or asymmetric. 



4 Hamel's magnetic top 

In this section an observed magnetism induced spinning phenomenon is analyzed. Up to 
our knowledge, the first magnetism induced spinning device was "Tesla's egg of colombus" 
exhibited in 1892 [11]. Using two phase AC energizing coils in quadrature, Tesla placed a 
copper plated ellipsoid on a wooden plate above a rotating magnetic field. The egg stood on 
its pointed end without cracking its shell and began to spin at high speed to the amazement 
of the scientists who witnessed the experiment. This effect was caused by induced eddy 
currents on the surface of the ellipsoid. A related magnetism induced spinning phenomenon 
is the Einstein-de Hass effect in which the rotation of an object is caused by a change in 
magnetization [9]. In this effect the magnetic field causes an alignment of electronic spins 
and their angular momenta is transferred to the atomistic lattice. We will show below 
through an idealized model based on magnetostatics that the spinning of the top in Hamel's 
device is due to surface friction. 

Mechanism behind curious rotation The following observation was made on a simple 
mechanical system consisting of a magnetized top and ring as shown in Fig. 1.1. When the 
ring is held above the top within a certain range of heights, and then tilted, one observes the 
top transitions from a state of no spin to a state of nonzero spin about its axis of symmetry. 

The system is modeled as two rigid bodies in magnetostatic interaction: a ball with a 
magnetic dipole aligned to one of its axes, and a fixed magnetized ring of radially aligned 
magnetic dipoles. The main tool used to analyze the observed curious rotation is a Lagrange- 
Dirichlet stability criterion [16]. It is shown that the fixed points of the system's governing 
equations correspond to the magnetic top being at rest with its axis of symmetry aligned 
with the local magnetic field and its translational position at a critical point of the magnetic 
potential energy. Stability of this point is determined by analyzing the nature of this critical 
point. If the attitude of the ring is normal to the surface, the magnetic potential energy 
is very nearly axisymmetric. A cross-sectional sketch of this potential energy is shown in 
Fig. 4.1. In this case there exists a ring of minima that are not individually stable. However, 
if the ring is tilted slightly, the potential has a unique local minimum opposite a saddle point 
as shown in the same figure. 

If the ring is tilted and the ball's initial position is unstable, the ball moves towards the 
stable fixed point which induces a resisting frictional force (see Fig. 4.2). The torque due to 
the sliding friction is in directions orthogonal to the moment arm q. However, the torque 
due to the magnetic field counters the torques about the axes perpendicular to ^3. This 
magnetic torque keeps ^3 aligned with the local magnetic field. Thus, the torque due to 
friction mainly causes a spin about the ^3 axis. 

Thus far, the ring has been kept fixed. If the ring is perturbed slightly the position of 
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the local minimum will change. Hence this system is unstable with respect to perturbations 
of the ring. This instability will be utilized to design a prototype stochastic mechanicial 
rectifier. 

No Till 

w 

/ 

ring of minima 

Figure 4.1: Cross-Section of Potential Energy Surface. Cross-sectional sketches of the magnetic 
potential energy with and without tilt. With no tilt in the ring, the potential energy surface is axisymmetric 
and has a ring of minima. If the ring is tilted slightly, a unique local minimum exists across from a saddle 
point. 



With Tilt 




global minimum 




Figure 4.2: Mechanism behind spin. A sketch of the ball, its axis of symmetry $3, the ambient magnetic 
field (blue arrows), the restoring magnetic torques, and frictional force. The ball tends to a position that 
minimizes its potential energy. If the ball is initially unstable, then as it moves towards the minimum a 
frictional force resists this motion. The frictional force causes torques about axes orthogonal to the moment 
arm q. However, the restoring magnetic torque counters the torques about the axes orthogonal to ^3 as 
shown in the sketch. This argument clarifies why the spin is primarily along ^3. 




Figure 4.3: Illustration of Magnetized Rigid Ball and Ring. The figure depicts a rigid ball on a 
flat surface with a dipolo fixed at the ball's centroid Oi,aii and in the direction ^3. We will assume that the 
center of mass of the ball is coincident with Of,aii- The ring is located at a height Hring above a reference 
point O on the surface, and its radius is r^ing. The magnetic dipoles are placed symmetrically around this 
ring with dipole moments in the radial direction. 



Hamel's magnetic top is modeled as an axisymmetric rigid ball of radius r and mass 
m. A magnetic dipole is attached to the center of mass, and in the direction of the axis 
of symmetry ^3 € M'^ as shown in Fig. 4.3. For simplicity, the surface friction is modelled 
using a sliding friction law. In what follows the continuous and discrete model are derived 
using the Hamilton-Pontryagin (HP) variational principle [3]. 

Magnetostatic Field Let B : K'^ — > be the ambient magnetostatic field. For simplic- 
ity, we assume the magnetostatic field is due to a magnetic ring consisting of N magnetic 
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dipoles equally spaced around a ring of radius rring centered at the point Oring and in the 
plane defined by the vector ^3 e M^. The point Oring is at a height H^ing above a refer- 
ence point O on the surface of the ball. To the reference point O, an orthonormal frame 
{ex,ey,ez) is attached. For i — I, ■■■,N, let di G M.^ denote the location of the ith-dipole 
with respect to Oring and nij e denote its dipole moment. The dipole moments are 
assumed to be aligned in the radial direction of the ring. Let be the permeability of free 
space. The magnetic field of the ith-dipole at a field point x e can be computed from 
the vector rj connecting the field point to the ith-dipole given by 

= X f^ring^i ^ring^z 

using the standard formula for the magnetostatic field due to a dipole: 

B(r„mO = (3(r,Tm,)r, - |lr,f m,) . (4.1) 

Using the principle of superposition, the magnetic field at a field point x G due to the 
N dipoles is determined as the vector sum of the magnetic fields due to each dipole: 

N 

B(x) = ^B(ri,m,). 




(a) axisymmctric (b) tilted (c) tilted 

^R = 0,eR = cl)R = Tr/4,0R = O = n/4, Or = Tr 

Figure 4.4: Illustration of Magnetized Ring and its Magnetic Field. The figures depict various 
orientations of a magnetized ring consisting of a discrete number of dipoles equally spaced around a circle 
in the plane defined by C3 which is also depicted. The magnetic fields lines in the plane y = are plotted. 
The angles (pn and On shown are the altitude and azimuth of the normal to the plane in which the ring is 
in. The height of the center of the ring, the radius of the ring, the number of dipoles, and the orientation 
of the ring are all parameters in this magnetostatic field model. 



Lagrangian of Magnetic Ball The configuration space of the system is Q = M"^ x 
S0(3) and its Lagrangian is denoted L : TQ ^ R. Let (x(t),x(t)) e x denote 
the translational position and velocity of the ball measured with respect to the reference 
frame attached to O. Let (ei, e2, e^) denote an incrtial orthonormal frame attached to Obaii 
and related to a body-fixed frame (^1,^2; €3) via the rotation matrix R{t) G S0(3). Let 
J = diag(/i, /2, /a) be the standard diagonal inertia matrix of the body and ^3 the axis 
of symmetry of the ball. We assume the mass distribution of the body is symmetric with 
respect to the axis of symmetry ^3 . This assumption implies that the principal moments of 
inertia about and ^2 are equal, i.e., I — Ii = l2- Let (x(t),x(t)) e x denote the 
translational position and velocity of the ball. 

We will use the isomorphism between and the Lie algebra of SO (3), so (3), given by the 
hat map (See appendix.). In terms of this identification, we define the reduced Lagrangian 
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^ : M3 X m3 X S0(3) X M3 ^ M as 

^(x,x,i?,a;) ^ L(x,x,i?,<I;i?). 
For the free magnetic ball, i.e., magnetic ball without dissipation, it is given explicitly by, 



^(x,x,i?,a;) = ^u;^ RJR^uj ~ mgx^es + (JB{x). 



(4.2) 



From left the terms represent the translational and rotational kinetic energy of the ball, the 
gravitational potential energy and the dipole potential energy. The ball is also subject to 
the surface constraint (/s : — > M given by: 



<^(x) = -r + GgX. (4.3) 
This holonomic constraint restricts the translational motion of the ball to a plane. 

Governing Conservative Equations The equations of motion will be determined using 
a HP description [3] . The constrained HP action integral is given by. 



dt. 



£{x, V, R, uj) + (p, X - v) + (^n, RR:^ + A(^(x) 

The HP principle states that 

Ss = 

where the variations are arbitrary except that the endpoints (x(a),i?(a)) and (x(6),i?(&)) 
are held fixed. The equations are given by. 



X = v (reconstruction equation) 
d£ 



dt' 



9v 

ox ox 



(Legendre transform), 
(Euler-Lagrange equations) , 



ip(x) — (constraint equation). 



dt 



R = uiR (reconstruction equation), 



dl 



(reduced Legendre transform). 



de 



R"^ — Tvu) (Lie-Poisson equations) . 



dt dR 

Evaluating these equations at £ as defined in (4.2) yields 



X = v 

mv = (DB(x))^|3 + (A-mg)e3 

T 

X e3 = r 

R = QR 

TT = RJR^ui 

iv = i^B{x). 



(4.4) 
(4.5) 

(4.6) 
(4.7) 
(4.8) 

(4.9) 

(4.10) 



(4.11) 
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Remark 4.1. Since the system is axisymmetric the Legendre transform in (4.11) simphfies: 

n = Iu:+ih- /)(u;T|3)e3 =^ = j (^n + ^-^{tt^ i,)^,^ . (4.12) 

As a consequence one does not need to solve for the evolution of all three columns of R{t) 
to integrate the ODE in tt. Instead one just needs to solve for the evolution of the third 
column, using: 

i3 = <2^3- (4.13) 
The following conservation law follows from axisymmetry. 

Proposition 4.1. The following momentum map is conserved under the flow of (4-11), 

Proof. This momentum map is due to an symmetry of the Lagrangian about the axis 
^3 which can be computed by formula (12.2.1) of [16]. The group acts on Q by: 

<i>?(x,i?) = (x,exp(s|3)i?). 

The corresponding infinitesimal generator is given by: 

V''3(x,i?) = ^ci>?(x,i?)|^^„. 

Thus, one can invoke Noether's theorem to conclude the following momentum map J : 
TQ -^{S^Y is preserved 



This conservation law indicates that if initially the top is not spinning about its axis 
of symmetry, then it will never spin about this axis. Thus, one cannot obtain the curious 
rotation using magnetic and gravitational effects alone. This result suggests surface friction 
will also play a role in producing this phenomenon. 



Governing Nonconservative Equations Let q = —re^ denote the vector connecting 
the center of mass C to the contact point Q as shown in Fig. 4.3. We model the surface 
frictional force using a sliding friction law proportional to the slip velocity, i.e., the velocity 
of the contact point on the rigid body relative to the center of mass Vg: 

F/ = -cVq 

This law assumes zero static friction but, is nevertheless reasonable on very slippery surfaces 
where static friction is negligible. Moreover, it is an experimental fact that the magnetic 
top exhibits this curious rotation even on oily surfaces. In reality one must keep in mind 
that the ball is not in point-contact with the surface; rather, a finite area of the ball is in 
contact with the surface and is moving relative to the surface to make spinning possible. In 
this case this sliding model of friction is quite reasonable. A more refined model of friction 
would include rotational torque and dry frictional effects (Coulomb friction). 
The slip velocity is given by: 

VQ=x + wq. (4.14) 
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(4.15) 



The force of friction is therefore, 

F/ - -cVq, 

and the torque due to friction is, 

Tf = qF/. 

The governing dynamical equations of the magnetic ball with friction are given by: 



TOV = (DB(x))'^ + (A- 7715)63 -c(v + re3a;), 

T 

^3 = '^^3' 

= 7^(-+^(-^l3K3), 

TV = ^3B(x) + cre^(v + r&3u;). 
Remark 4.2. The fixed points of (4.15) satisfy: 
x = 0, uj = 

TT = =^ ^3 = K2B(x), K,2 G M, constant, 
v = =^ (DB(x))^B(x) =0. 

The Lagrange-Dirichlet criterion can be used to analyze the stability of these states since 
the kinetic energy vanishes at these fixed points and since the dissipation is proportional 
to velocity [16]. In particular, one can invoke a classical theorem due to Thomson- Tait- 
Chetaev, to conclude that if the fixed point is potentially stable or unstable, then it remains 
stable or unstable after introducing arbitrary dissipative forces proportional to velocities 
[17]. By this criterion, the fixed points are stable provided that the equilibrium is a strict 
local minimum of the potential energy. The potential energy evaluated at this equilibrium 
point is given by Vg : M'^ ^ M, 

14(x) = -K2B(x)TB(x). 

In the limit as the number of dipoles along the ring is infinite, B is axisymmetric with 
respect to the surface provided the attitude of the ring is normal to the surface on which 
the ball is on. In this case Ve can be written as a function of the distance to the origin: 
f{p) = Ve(x) where p = |jx|j with x"^e3 = r. Assume that the ring is above the ball and the 
potential energy is Ve(x) < for all x e M^. The origin is an unstable critical point, since 
/ has a local maximum at that point. Moreover, /(p) tends to zero as p becomes large. 
One can pick a distance M 3> rrmg (the radius of the ring) sufficiently large to make f{p) 
arbitrarily close to zero. In the interval p E {0,M), Ve{p) is continuous, and therefore, the 
function has a local minimum in this interval. However, since / is axisymmetric this critical 
point is not a local minimum in the plane, but rather a circle of critical points. As shown 
in Fig. 4.5, the ball will be unstable with respect to circumferential perturbations. The 
figure also shows that when the axisymmetry of the magnetic field is broken by changing 
the orientation of the ring, there is a critical point that is a local minimum of Ve(x). 



Simulations Four simulations are performed on the top to confirm the theory and explain 
the curious rotation of the physical system. We use the HP integrator introduced in [3]. Let 
ipB and 9b denote the azimuth and altitude of the ball's attitude ^3. In all of the simulations 
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(a) axisymmetric (b) tilted (c) tilted 

= 0, 0fl = 4>B. = ^/16, Or = 4>R. = ^/16, Or = TT 



Figure 4.5: Surface plot of Ve- The potential energy surface is plotted as a function of x and y for ball 
radius r = 0.3, ring height H = 2, 2b dipoles on the ring, and ring radius r^ing = 4. The angles <f) and 9 
shown are the altitude and azimuth of the normal to the plane in which the ring is in. By tilting the ring 
one can obtain a critical point which is a local minimum of Ve- 



the top is initialized as follows: 

c= 0.1 kg • m/s,r 1.8 cm,™ = 200 g,0B = 0,6*5 = 0,/3 ^ h = 2/5mr^ 
/x||m,||/(47r) = Ai||^3ll/(47r) = 10"^ Tesla meter^, 
Hring — 20 CHI, = 20 dipoles, rring — 34 cm 

The timestep size and number of timesteps are h = 0.025 and N = 20000 respectively. In 
the first four simulations the orientation of the ring is kept fixed. The first case simulated 
is when the magnetostatic field is axisymmetric and friction is absent. In this case the top 
undergoes what appears to be chaotic motion involving a balance between the magnetic 
potential and translational kinetic energies as shown in the simulation and Fig. 4.6. In the 
second simulation the altitude of the attitude of the ring is slightly perturbed causing the 
magnetostatic field to be noticeably asymmetric, see Fig. 4.7. However, in both of these 
cases the top does not acquire spin about the symmetry axis ^3 as predicted by the theory. 

In the presence of friction c = 0.3, the motion changes. The third case considered is 
the same configuration as the first case, but with friction. In this case the kinetic energy 
of the top is dissipated and the top moves towards a point where the magnetic potential 
energy is minimum, see Fig. 4.8. However, no spin develops. If the attitude of the ring is 
slightly perturbed, spin does develop as shown in Fig. 4.9. The origin of this spin in case 4 
is explained here. 

It is very clear from (4.15) that one can get spin about the axis of symmetry in the 
presence of friction. However, the frictional forces may produce torques about the other axes 
as well. So what is not clear is why the ball spins primarily about ^3. This phenomenon is 
clarified in the following remark. 

Remark 4.3. If the initial position of the ball is not a local minimum of the magnetic 
potential energy (as in the simulation), the position will be unstable as predicted by a 
Lagrange-Dirichlet criterion. The ball's position then adjusts to minimize its magnetic 
potential energy which causes sliding friction (see Fig. 4.2). The torque due to the sliding 
friction will introduce a torque in directions orthogonal to the moment arm q. However, the 
torque due to the magnetic field will counter the torques about the axes perpendicular to 
^3. Keep in mind that the magnetic torque keeps ^3 aligned with the local magnetic field. 
Thus, the torque due to friction mainly causes a spin about the ^3 axis as shown in the 
simulation. Snapshots of the simulation are provided in Fig. 1.2. 
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(a) potential surface Ve 



(b) momentum map J(t) 



Fi gure 4.61 Conservative Axisymmetric. The figures plot x{t\y{€) and J{t) of the top in the absence 
of friction, and for (pn = and Oji = 0. Superimposed on (a) is the potential energy level set which appears 
axisymmetric. (a) shows that the top bounces between the local maximum at the center and the circle of 
minima (dark blue ring), (b) shows that J{t) is preserved as predicted by the theory. The accompanying 
simulation vividly illustrates this motion and demonstrates that one does not get curious rotation in this 




(a) potential surface Ve 



(b) momentum map J{t) 



Figure 4.7: Conservative Tilted. The figures plot x{t),y(t) and J{t) of the top in the absence of 
friction, and for (f>ji = and Oji = 7r/64. Superimposed on (a) is the potential energy level set which is no 
longer axisymmetric. (a) shows that the top moves erratically within an annulus. (b) confirms that J{t) 
is preserved. The accompanying simulation vividly illustrates this motion and demonstrates that one does 
not get curious rotation in this case. 




(a) potential surface Ve 



(b) momentum map J{t) 



Figure 4.8: Nonconservative Axisymmetric. The figures plot x{t),y(t) and J(t) of the top in the 

presence of friction, and for iji^ = and 6^ = 0. Superimposed on (a) is the potential energy level set 
which appears axisymmetric. (a) shows that the top moves from the local maximum at the center to a 
point that minimizes the potential energy (dark blue ring). However, (b) shows that no spin develops. The 
accompanying simulation vividly illustrates this motion and demonstrates that one does not get curious 
rotation in this case. 
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(a) potential surface (b) momentum map J(i) 

Figure 4.9: Nonconservative Tilted. The figures plot x(t),y{t) and J{t) of the top in the presence of 
friction, and for = and Sjj = 7r/64. Superimposed on (a) is the potential energy level set which is no 
longer axisymmetric. (a) shows that the top moves along arcs, (b) confirms that one does get spin in this 
case. The accompanying simulation vividly illustrates this motion and clearly demonstrates that the ball 
goes from a state of rest to a state of nonzero spin about its attitude. 



5 Fluctuation Driven Magnetic Motor 

The instability of the top with respect to perturbations of the ring is the key idea behind 
the fluctuation driven magnetic motor. These random fluctuations are modelled as a white 
noise torque on the ring. However, if the ring is allowed to move freely, it could possibly 
turn on its side. To stabilize the ring, a fixed magnetized outer ring is installed. 

By adjusting the radii and heights of the rings and the inertia of the inner ring, one 
can obtain a configuration in which the attitude of the inner ring can be randomly torqued 
without undergoing large excursions from the vertical position. In this case numerical exper- 
iments reveal that one can adjust the amplitude of the white noise so that the ball undergoes 
directed motion on certain time-scales. To be precise the numerics indicates that starting 
from a position of rest initially the top's motion is dominated by the effect of white noise. If 
the outer ring is close enough one can see the top oscillate between the wells in the magnetic 
potential caused by the dipoles in the outer ring. This behavior is reminiscent of stochastic 
resonance. 

After some time, the top accumulates enough kinetic energy that it displays directed 
motion along a circle of certain radius. This motion is inertia-driven and the energy injected 
into the system by the white noise mainly adds to the speed of the top. Provided that the 
inner ring does not turn on its side, one of two things can happen: 1) the top reaches 
a critical velocity in which the amount of energy dissipated by the surface friction is on 
average equal to the amount of energy injected by the thermal noise or 2) the top gathers 
enough kinetic energy to escape from the potential well created by the inner and outer rings. 
Conducting this same experiment at uniform temperature, reveals that this phenomenon 
persists. The isothermal, magnetic ball-ring system is a prototype fluctuation-driven motor. 

The fluctuation-driven motor considered in this paper is related to the granular, mag- 
netic balls in ferrofluidic thermal ratchets. In such ratchets a time-varying magnetostatic 
field transfers angular momentum to magnetic spherical grains in a ferrofluid [8] . The ferro- 
magnetic grains are modelled in the same way as the rigid ball in the prototype. However, 
our work is different in that the magnetic potential in the prototype is autonomous. In 
fact, the use of a non-autonomous potential to design a fluctuation driven motor is well 
understood [1]. 

In this section the equations of motion for the dynamics of a magnetic ball interacting 
with a dynamic, inner magnetized ring are derived. The magnetic effects of a fixed outer ring 
are also taken into account. As before the magnetic ball is assumed to be constrained to a 
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flat surface that resists the motion of the top via surface sHding friction. The center of mass 
of the inner ring is kept at a fixed height, but otherwise it is free to rotate and subjected 
to white noise torques. The governing stochastic differential equations are analyzed using 
energy arguments. Simulations validate this analysis and show that one can get sustained 
directed motion of the top from random perturbations of the inner ring. Moreover, it is 
shown that this phenomenon persists even when the system is at uniform temperature. 
This system is called a fluctuation driven magnetic motor. 

The mass and radius of the magnetic ball are m and r respectively. The mass of the 
inner ring is mi„nor- The heights and radii of the inner and outer rings are plotted and 
defined in Fig. 5.1. 



Q,.„,- 





















c 





Figure 5.1: Illustration of Magnetized Rings and Magnetic Top. The figure shows the magnetic 
top on the surface with two rings above it. The inner and outer rings have radii router and Hnner S'ld 
heights -Houtor and Hinncr respectively. The centroids of the ball, inner and outer rings are also labeled. 
The purpose of the outer ring is to prevent the trivial equilibrium in which the inner ring turns on its side. 



Lagrangian of Magnetic Ball &: Inner Ring The configuration space of the system 
is Q = M'' X S0(3) X S0(3). Let Obaii, Oinner and Oouter denote the centroids of the 
ball, inner and outer rings respectively. Let (ei,e2,e3), (fi,f2,f3), and (gi,g2,g3) denote 
inertial orthonormal frames attached to Obaih Oinner and Oouter, and related to body-fixed 
frames (^1,^27^3)- (Ci,C2.C3)- and {Vi,V2,V3) via the rotation matrices Rsit), Rnit), 
and Ro G S0(3). Let 1b — diag(/i, /2, 13) and 1;^ = diag( Ji, J2, J3) be the standard 
diagonal inertia matrices of the body and inner ring, and let ^3 and ^3 be the attitude 
of the ball and inner ring. The following assumptions are made: the outer ring is fixed, 
the mass distribution of the inner ring is symmetric with respect to its attitude (or axis 
of symmetry), and the ball's mass distribution is spherically symmetric. This assumption 
implies that Jb = h — I2 — I3 and J = Ji = J2. 

The magnetostatic field at any field point is due to the dipole in the ball, and an inner 
and outer ring consisting of N and M magnetic dipoles respectively. On each ring the 
dipoles are equally spaced around a circle of radius Tinner and roster in the planes defined 
by the vectors 773 and ^3 e M.^ . We assume that there is no self-interaction between the 
dipoles within each body. The heights of the rings above the surface are denoted by Hi^ncv 
and -ffoutcr- 

For i = 1, • • • , iV and j = 1, • • • , Af, let "'^'■(i) and d°"'"' e denote the location of 
the ith and jth-dipoles on the inner and outer rings with respect to the points Oinner and 
Oouter and let mj""°''(t) and m°^''''' G M.^ denote the orientation of their respective dipole 
moments. The dipole moments are assumed to be in the radial direction of the ring. The 
magnetic field of each dipole at a field point can be determined from the vector r connecting 
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the field point to the dipole using (4.1) 

B^""°'(r) = B(r, mj""") (field of inner ring dipole), 
gouter^^) ^ m°"*") (field of outer ring dipole), 
Bo(r) B(r, ^3) (field of ball). 

Define the following vectors Vij G M.^ i — 0, - ■ ■ ,N and j = 0, - ■ ■ ,M which joins the inner 
and outer dipoles, the inner dipoles to the ball, and outer dipoles to the ball as follows, 

_ joutcr „ Hontcres + dj""''^ + 77i„„ere3 if 1, j > (ith outer to jth inner dipole) , 
X - dj""" - HinncrGs if j = (jth outcr to ball) , 

X - d°"*"' - iJoutcrCs if i (ith inner to ball). 

The spatial representation of the reduced Lagrangian i : TR^ x S0(3) x E-^ x S0(3) x R"^ 
M of the free magnetic ball, i.e., magnetic ball without dissipation, is given by, 

N M N M 

+ 2 ^ € jBr-(r.o) + E ljB°"*-(ro,) + ^ ^(mr-)TB°"*- (r,,) (5.1) 

i=l 3 = 1 i=l j = l 

From left the terms represent the translational and rotational kinetic energy of the ball, the 
rotational kinetic energy of the ring, the gravitational potential energy of the ball and the 
magnetic potential energy of the inner ring and ball dipoles. 

Governing Equations for Magnetic Ball &; Ring The equations of motion are de- 
termined from a Hamilton-Pontryagin (HP) principle [3]. The HP action integral is given 
by, 

s = J ei-x,±,RB,u;B,RR,i^R) + (t^b,RbR^ - i^") + (t^r,RrRr dt. 
The HP principle states that 

Ss = 

where the variations are arbitrary except that the endpoints {x{a), RB{a), RR{a)) and 
{x{b) , RB{b) , RR{b)) are held fixed. The equations are given by, 

d de _ d£ ^dif 
dt dx dx dx 
ip{x) = 

— i?B = i^bRb 

dt 

-j^Rr = ^rRr 

de 



dujB 

de 



d ^ de ^=r~~ 
tt^b 



dt dRB 
d ^ _ di 
Jt'^^ ~ Mr 



(Euler-Lagrange equations for ball). 


(5.2) 


(constraint equation) , 


(5.3) 


(reconstruction equation for ball). 


(5.4) 


(reconstruction equation for ring). 


(5.5) 


(reduced Legendre transform for ball). 


(5.6) 


(reduced Legendre transform for ring). 


(5.7) 


(Lie-Poisson equations for ball). 


(5.8) 


(Lie-Poisson equations for ring). 


(5.9) 
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The nonconservative system is obtained by adding the frictional force and torque derived 
earher to the translational and rotational equations of the ball: (5.2) and (5.8). 

Random Perturbations &; Uniform Temperature Consider driving the inner ring by 
the following Wiener process. Let W e M'^ denote Brownian motion in M.^, and append the 
following random torque: 

Tr ^ adW e 

to the Lie-Poisson equation for the ring (5.9). 

The work performed by T/j is equal to the change in total energy of the system. On the 
other hand, the work done by the kicks over a time interval [a, b] is given by the following 
formula: 

Work of kicks = a / T^wj^, (5.10) 

J a 

The work transferred from the inner ring to the top is given by: 

Work transferred to top = ATotal Energy — ARing Energy. (5.11) 
Thus a measure of the efficiency of the magnetic motor is given by the following ratio: 

„ Work transferred to top , , 

- Work of kicks ■ 

As a next step frictional torque is introduced to the inner ring and thermal torque (white 
noise) to the ball, so that the generator of the process is characterized by a unique Gibbs- 
Boltzmann invariant distribution. The reader is referred to the Appendix for the governing 
equations of the magnetic motor at uniform and non-uniform temperatures. 

Simulations A stochastic variational integrator is used to carry out these simulations [5]. 
Simulations are conducted at uniform and non-uniform temperature as described below. 
The initial conditions and parameters used are given by 

c = 0.15 kg • m/s, r = 4.0 cm, m = 500 g, (j)B = 0,03= 0, /g = /i = 2/5mr^, 

^^ll^inncr, out.r y ^^^^ ^ fMsW/i^^) = 2 X IQ-' T • m^ 
-ffinncr = -Pouter = 48 Cm, = 20 dipolcS, M = 5 dipolcS, 
nnnor = 75 Cm, router = 98 Cm, TOi„ncr = 6 kg, J3 = minncrf fnncr ' J = 1/2 Js- 

In all of the simulations the magnetic top is initially at rest with its attitude aligned to 
the vertical. The difference in the simulations is the amplitude of the oscillations and the 
discrete sample from the normal distribution. Figures 5.6-2.2 present data from the uniform 
temperature simulations. The key point about the uniform temperature simulations is as 
seen in Figures 5.9 and 2.2 the magnetic top transitions from noise to ballistic motion. As 
expected the non-uniform temperature case also exhibits this behavior as shown in Figures 
5.2-2.1. 
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Figure 5.2: Total and inner ring energy (non-uniform temperature). The total energy of the 
system is shown in green and the ring energy in cyan for different amplitudes of the noise and for two 
different samples. Although there is dissipation, the thermal fluctuations inject energy into the system. 




(a) a = 0.0002 (b) a = 0.0002 (c) a = 0.00025 

Figure 5.3: Energy injected and dissipated (non-uniform temperature). A plot of the energy 
injected by the white noise as computed analytically and dissipated by friction computed numerically. 




(a) a = 0.0002 (b) a = 0.0002 (c) a = 0.00025 

Figure 5.4: xy-position of magnetic ball (non-uniform temperature). An aerial view of the path 
of the center of mass of the ball in the xy-plane. 




(a) a = 0.0002 (b) a = 0.0002 (c) a = 0.00025 

Figure 5.5: xy-position of magnetic ball planar view (non-uniform temperature). The x and y 

components of the center of mass are plotted in green and blue respectively, ft is very clear from this plot 
that the motion transitions from noise driven to being inertia driven. 
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(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 

Figure 5.6: Total and inner ring energy (uniform temperature). The total energy of the system 
is shown in green and the ring energy in cyan for different amplitudes of the noise. Although there is 
dissipation, the thermal fluctuations inject energy into the system. 




(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 

Figure 5.7: Energy injected and dissipated (uniform temperature). A plot of the energy injected 
by the white noise as computed analytically and dissipated by friction computed numerically. In the constant 
temperature case we see an approximate balance between these energies. 




(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 

Figure 5.8: xy-position of magnetic ball (uniform temperature). An aerial view of the path of the 
center of mass of the ball in the xy-plane. 




(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 

Figure 5.9: xy-position of magnetic ball planar view (uniform temperature). The x and y 

components of the center of mass are plotted in green and blue respectively. It is very clear from this plot 
that the motion transitions from noise driven to being inertia driven. 
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6 Conclusion 

This paper introduces a novel mechanism to rectify random perturbations to achieve directed 
motion as a mcta-stablc state and ballistic mean-squared displacement with respect to the 
invariant law. The basic idea behind the mechanism is explained using the simple sliding 
disk model. With this model one can prove that if the potential energy is non-constant, 
then the invariant Gibbs measure of the system, /i, is ergodic and strongly mixing. As 
a consequence /x a.s. it is shown that the translational displacement of the sliding disk is 
not ballistic. However, it is shown that the mean-squared displacement with respect to the 
invariant law is ballistic. 

This anomalous diffusion manifests itself numerically. That is, starting from a position 
at rest and when U is non-constant, one observes flights in the ai-displacement over time- 
scales which are random and the probability distribution of the flight-durations have a heavy 
tail. Moreover, from a numerical statistical analysis, it is observed that the mean-squared 
translational displacement of the sliding disk is ballistic and its translational displacement 
is characterized by long time memory /correlation. The paper proceeds to show that this 
basic phenomenon arises in a more complex, rigid-body system consisting of two rigid bodies 
interacting via magnetostatic effects. Along the way we explain the observed dynamics of 
Hamel's magnetic device. 

7 Appendix 

so (3) Preliminaries In the paper we use the following isomorphism between an element 
of the Lie algebra of S0(3), Tg SO (3) = so (3), andM"^. Recall that elements of so (3) are skew- 
symmetric matrices with Lie bracket given by the matrix commutator. Let w = (wi, W2, t^s)- 
Then one can relate with a skew-symmetric matrix via the hat map ^: — > so(3). 








-UJ3 


UJ2 


Q = 


0.-3 





-UJl 




—UJ2 


UJl 






Let g{t) be a curve in S0(3). With this identiflcation of so(3) to R^, the right-trivialization 
of a tangent vector g to this curve, given by ^ = TRg-i • 5 e so(3), can be written in terms 
of the spatial angular velocity vector u> e M^, i.e., as ^ = Q G so(3). 

Magnetic motor at non-uniform temperature governing equations Below the gov- 
erning equations of the ring-ball system with nonconservative effects due to surface fric- 
tion and white noise are written in Ito form. We introduce the potential energy function 
[/ : X S0(3) X S0(3) M which represents the total potential energy of the ring-ball 
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system. In terms of U the governing equations can be written as, 



dx 


=vdt 


dv 


1 f dU 
m \ 9x 


T 
X 63 


— V 




=ujBRBdt 


cIRr 


=ijJJiRRdt 




=Jb<^b 




—RR^RRjii^R 


dins 


"I dRs"^^ 


d-rf^ 


V ORr 



+ (A — mg)e^ — cVq j dt 



Rb - cqVg dt 



Evaluating these governing equations at ^ defined in (5.1) yields 
dx =wdt 

. ( N M \ 

^^^^ I 2^DBr-(r,o)^43 + E'^B-*-(ro,)T|3 j dt 

+ - (-cVq + (A - mg)e3) dt 
m 

T 

X 63 = r 
dRs =CjBRBdt 
dRR =CJ~RRRdt 

T^B — JbI-^B 

TTfl —RrIrR^Ur 

( N M \ 

dTTB =f3 2 E Br"-(r,o) + B°"*"(ro,) dt - cqVgdt, 



i=l 



3 = i 



N 



dt 



diTR =2^ [H4^^Bo(r,o) - (DBr-(r,o)^€3) 

N M ^ ^ 

,innor-iDOutcr /„ \ , jinncr /'■r»T50utcr/„ ^T^inncr^ 



- E [mr""^B°-- (r,,) + dj-- (DB°--(r,,)-m5' 

i=i j=i 

adW. 



dt 



(7.1) 



(7.2) 

(7.3) 
(7.4) 
(7.5) 
(7.6) 
(7.7) 

(7.8) 



(7.9) 



Since the ring is axisymmetric its Legendre transform simplifies: 



J- J: 



(7.10) 



TVR = Jur + ( J3 - J)(a;^C3)C3 =^ i^R= J [t^r 
Remark 7.1. The energy of the magnetic motor is given by 

, . in rp 1 IT "^T IT / \ 

-E(x,V,i?B,7rB,i?ii,7r_R) = yV V+-J^ TVgTTB + -TTrRrIr RRTTR + Uili.RByRR)- 
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Using Ito's formula one can calculate the stochastic differential of the energy in the case 
when the magnetic motor is at non-uniform temperature 

dE = ^-cV^Vq + trace [l^^] ^dt + aw^tfW. 

Integrating yields, 

/■* 1 

E{t) - £'(0) = -c / V^Vgds + -a^ trace [l^^] t + martingales. 

The first term represents the energy dissipated by friction. The next terms represent the 
work done by the white noise torque as computed by Ito's integral. 



Isothermal, magnetic motor governing equations Similar to the sliding disk, to put 
the magnetic motor at constant temperature we define the following dissipation matrix. 



C = 



l/m^ 

l/m^ r/(mJ) 

r/(mJ) rVj2 
-r/(mJ) 



— r/(mJ) 





The translational position of the ball is written in coordinates as x = {xi,X2,r). The 
translational and angular velocities of the ball are given by v = (fi,i'2,0) and ojb = 
(wg^ w^^o;^^). The dynamical equations for the constant temperature magnetic motor 
are given by: 



dxi 

dx2 

dRg 

dRR 



v\dt 

V2dt 

ajBRsdi 
uTkRiidt 



d,vi 




Ujei/m 




' mvi ' 




dBy-^ 


di'2 




f/Jea/m 


dt — cbC 


mv2 


+ asC ' 


dBy^ 
dB 


du;'-^\ 




U^e2/J_ 








dB (2) 



dto 



(3) _ 



B 

dnji 



Rr^rRJ{<^r. 

~Ur - CrU)r 



ctRdBR 

The terms Ux, Ub, and Ur are defined in terms of the inner product on 

'dU 



as, 



dx ' ' 
dU 

dU_ 

R 



dR 



Rl,y 



Rl,y 



(7.11) 

(7.12) 
(7.13) 
(7.14) 

(7.15) 

(7.16) 
(7.17) 
(7.18) 



and dR^U,dR^U : S0(3) ^ T* S0(3), and d^U : I 
uniform temperature, one can prove the following. 



dxU{x,RB,RR) ■ y 

Qrb U {x, Rb, Rr) ■ yRs 
dRjJ{x,RB-,RR) ■ yRR 

' T*m?. Similar to the sUding disk at 
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Theorem 7.1. Let x denote the phase space of the magnetic motor. Set ^ to be the solution 
of (7.11)-(7.18). Suppose that Cfi/aj^ = cb/o:%- Let E denote the energy of the magnetic 
motor given by, 

£;(x,v, i?B,-n-s,i?i^,7r7j) = — V v+-J^ tt^ttb + -Tr^ji^fll^ Rj^tvr + U{x, Rb, Rr) ■ 
Let [3 — 2cb/c(%- The following measure, 




Z 



(7.19) 



is the unique Gibbs invariant measure of the stochastic process ^. 



Remark 7.2. By using Ito's formula one can show that: 



dE 



( 



( 



m 



1 



H — j\ dt + -afj. trace [ij^./^] ) dt + martingales. 



Integrating yields, 



E{t) - E{0) 
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